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Abstract 

In the case of an external Hamiltonian at most quadratic in position 
and momentum operators, we use the ABCD^ formulation of atom optics 
to establish an exact analytical phase shift expression for atom interfer- 
ometers with arbitrary spatial or temporal beam splitter configurations. 
This result is expressed in terms of coordinates and momenta of the wave 
packet centers at the interaction vertices only. 



1 Introduction 



Recently atom interferometers [Q have been described by the ABCD^ formahsm 
of Gaussian atom optics , which yields an exact formulation of phase shifts 
taking into account the wave packet structure of atom waves. 

For the theory of atom interferometers two basic stages are required: 

1. a proper description of the propagation of wave packets between the beam 
spHtters 

2. an adequate modelization of the beam sphtters themselves. 

The first stage is achieved through the ABCD^ theorem whose main results 
are briefly recalled in section 2. The second problem is adressed by the ttt 
theorem which provides a simple model for the phase introduced by the splitting 
process. 

In this paper we give a compact way to express the atom interferometer phase 
shifts in terms of the coordinates and momenta of the wave packet centers only. 
For this purpose we derive two new theorems (the four end-points theorem and 
the phase shift theorem) valid for a Hamiltonian at-most-quadratic in position 
and momentum operators. 
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2 The ABCD^ theorem 



In this framework we consider a Hamiltonian which is the sum of an internal 
Hamiltonian Hq (with eigenvalues written with rest masses rrii) and an external 
Hamiltonian Hext'- 

Hext = '^Pop^ (*) ■P^~-^<l^-7 (t) -(hp- ^ (t) ■ i<h^ X P^p)-'m~g (t) .q^ (1) 

where one recognizes several usual gravito-inertial effects: rotation through 
(t), gravity through ~g (t), gradient of gravity through 7 (t),... and where 
g (t) is usually equal to the unity tensor (for simplicity wc omit the transposition 
sign ~ on vectors). 

For a wave packet tp {q,ti) — wp {q,ti, qi,pi, Xi,Yi), where qi is the initial 
mean position of the wave-packet, pi its initial mean momentum, and (Xi, Yi) 
its initial complex width parameters in phase space, one obtains the ABCD^ 
theorem 

i/j{q,t2) = J d^q'.K{q,t2,q\ti).wp{q',ti,qi,pi,Xi,Yi) 

= e*''^'^''^''^'''^'"'>.wp{q,t2,q2,P2,X2,Y2) (2) 

where K and Sd are the quantum propagator and classical action respectively, 
and where q2,P2, X2,Y2 obey the ABCD law (G and U are the representative 
matrices of (t) and the rotation fl (t), and we write A21 instead of A (<2,ii) 
for simplicity): 
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For example the phase of a gaussian wave packet is: 
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Scl{t2M,qi,Pl) lf^+P2-{q - q2) /n+^iq - q2) .'Rc{Y2.X^ ^).{q-q2) (5) 
and in this case the main phase shift recorded between ti and ^2 is equal to: 

Sci {t2,ti,qi,pi) /h+pi.qi/h- P2-q2/TT' (6) 
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3 The ttt theorem 



When the dispersive nature of a laser beam splitter is neglected (wave packets 
structure preserved), its effect may be summarized by the introduction of both 
a phase and an amplitude factor (see jl^ and ji) for a detailed proof): 

Mf,a.e-*('^**-'^*-'*+'^') (7) 

where t* and q* depend on Ia and qa, the mean time and position of the laser 
wave. 

For a temporal beam splitter: 

t* = tA 

q* = Qci (tA) 
k* = k 

ip* = If (laser phase) (8) 
For a spatial beam splitter: 

q = qA 

t* such that qd (i*) = qA 

k* = k + Sk 

ip* = ip + Sip (9) 

where 5k is the additional momentum transferred to the excited atoms out of 
resonance, and where 6ip is a laser phase: 6ip = —Sk.qA (see [||). 

Let us emphasize that these calculations do not rely on the assumption that 
the splitter is infinitely thin or that the atom trajectories are classical. 



4 The four end-points theorem for a Hamilto- 
nian at most quadratic in position and mo- 
mentum operators 

We shall cut any interferometer into as many slices as there are interactions on 
either arm and thus obtain several path pieces (see section 5). From now on we 
shall consider systematically pairs of these homologous paths (see Fig. 1) in the 
case of a Hamiltonian at most quadratic. 
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Fig. 1: A pair of homologous paths 

These two classical trajectories are labelled by their corresponding mass {ma 
and mp), their initial position and momentum {qai, Pai, 9/31 and ppi) and their 
common drift time T — t2 — ti. 

Before establishing the first new theorem let us consider the expression of 
the classical action for the a path (see 

Scl{t2,h,qal,Pal) cTj n-'^ I A i D / \ : ^ ^-^ 

= ^.U.G .[A.qai+B.paima)+ at 

via Jti ™a 

AC Pal BD Pal Pal ,-,„s 

+qai-—;r-qai H —^r- 1 -BC.qal (10) 

2 nia 2 rua rua 

which can be rewritten as: 

Scl{t2,tl,qal,Pal) Pa2 Pal ^ i: r<-l TT C , f^' 

= 7, ■<lo^2--7, -qal ~ 7^£.-U .G .U.^+ dt 

rUa ZrUa Inia 2 Jf^ ma 

+ hu.G-\qa2-U.U.^-^ (11) 

z 2 TOq 

with the help of the definition of qa2 and (see ||). Then we can use the (i 
path to replace S^.U.G^^ with: 

lU.G-^^^-^-C.qpi^D.P-^ (12) 

mf3 TO/3 

Consequently we get: 

Scl{t2,tl,qal,Pal) 1 fPa2 , P/32 \ 1 f Pal , P/31 \ 

^ ■qa2 - 77 1 -qal 



ma 2 \TOq, To/3 / 2 \ma TO/3^ 

+/^(t2,<i) +/(a,/3) (13) 

where h (^2, ^i) is independent of positions and momenta and where / (a, /3) = 
f {f3,a). The same goes for the expression of Sd it2,ti,qfii,Pi3i) /"^/3 which is 
obtained by exchanging a and /?. Finally we arrive at the first new theorem (a 
more general demonstration starting with Hamilton principal functions is given 
in appendix 1): 
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Theorem 1 

Sclit2,tl,qal,Pal) ^ f Pa2 , P/32 \ 
_ Scl{t2,ti,qpi,Pf3i) I fpa2 
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or equivalently: 
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(15) 



which gives the main part of the phase shift expressed with the half sums of the 
coordinates and the momenta of the four end-points only. 

In the case of identical masses {nia = mp) this expression simplifies to: 



Scl {t2,tl,qal,Pal) -Pa2-qa2 +Pal-9al 
- [Scl {t2,tl,qpi,ppi) -p/32.9/32 + P/31 -9/31] 

'9a2+9/32\ / N (qal+qpi 



(P/32 -Pa2) 



(P/31 - Pal) • 



(16) 



5 The phase shift theorem for a Hamiltonian 
at most quadratic in position and momentum 
operators 

In this section we draw on the results of previous sections to establish the inter- 
ferometer phase shift expression for an arbitrary beam splitters configuration. 

For a sequence of pairs of homologous paths (an interferometer geometry) 
(see Fig. 2) 
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Fig. 2: Interferometer geometry sliced into pairs of homologous paths between 
interactions on either arm (when an interaction occurs only on one arm the 
corresponding k on the other arm is set = 0) 

one can infer the general sum for the main coordinate dependant part of the 
global phase shift: 



Pl3D - PaD ( ^ _ qaP + qi3D \ Pal + P/31 



N 



2h 



qai + qpi 
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(17) 

If now wc take into account the other terms of the phase shift we finally 
obtain the following result (given here for a Gaussian wave packet) : 

Theorem 2 

^(l){q,tN+l=tD) = {P0D - Pao) ■ {q - ^"^ ^ ^ 1^- 2^^^^ ' ^^^^ ~ ^"^^ 

^ r ■+ ■ 



i=l 
N 



+ 



~ + ~ ^"•'+^^ 2n • ^^'^^ ~ 

2% 

ma,N 

2% 



{q - qpo) ■ Re {Yd.Xjj^) . {q - quo) 
{q - qao) ■ Re (Yd.Xj^^) . {q - q^o) 



(18) 



where S^i = {ti+i,ti,qai,Pai + hkai,mai). 
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This fundamental result is valid for a time-dependent Hamiltonian and takes 
into account all the mass differences which may occur. It allows to calculate ex- 
actly the phase shift for all the interferometer geometries which can be sliced as 
above: symmetrical Ramsey-Borde (Mach-Zehnder), atomic fountain clocks,... 
All these particular cases will be detailed in a forthcoming paper (see Q) 

Let us point out that the nature (temporal or spatial) of beam splitters 
leads to different slicing of the paths. In the spatial case, indeed, the number 
of different t* may be twice as great as in the temporal case (see the definition 
of t* in these two different cases in section 3). 



6 Phase shift after spatial integration 



In an actual interferometer one has to integrate spatially the output wave packet 
over the detection region. With Gaussian wave packets this integration leads to 
a mid-point theorem |^ " The first term of Acj){q^t£,) disappears when the 
spatial integration is performed" . 

Furthermore the terms which depend on the wave packets structure {Y and 
X) vanish when mp,^ — ma^N (which is always the case). One obtains finally: 



A(f> {to) 



Pal + Ppi 



2h 
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{( Z— + TTZ — ■ - QaA+l) • {q|3^ - qai) 



S^i Pl3,i+1 , ^ 

r T; • Vla.i+l — qf3,i+l) 

171/31 2mf3i 
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7 Identical masses and symmetrical case 

The case of identical masses is an important approximation which is commonly 
used for the modelization of many devices like gravimeters and gyrometers W , 
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[§i @- If TUai — Tupi — TTi, Vi, this general phase shift becomes : 

l\4>[tD) = 7^ • (g/31 - gal) + 2^ (K/3i - fcai) ■ ^ 

1=1 

i=l 

We can also specify the form of this phase shift when the interferometer 
geometry is symmetrical (see Fig. 3). 



Fig. 3: A typical symmetrical interferometer 

This symmetry is expressed as : kfn + k^i = 0, Vi G [2, iV — 1], i.e. it is 
a symmetry with respect to the direction of the particular vector : Pinuiai + 

^kijiitial/'^- 

Consequently : 

N-l N-1 

A4> (tN) = fci.gi + 2 y: kj-^^ + k,.^-^i^L±^ - y: {^,. - (21) 

i=2 i=l 

But Vi e [2, TV - 1] : 

= Wl.l + ^J+l,!--?! H -(^1 + ^) 

m 2 

= Q(t.+i) (22) 

which can be calculated with the ABCD^ law. 

It depends only on qi (initial position) and pi + ( "Bragg initial momen- 
tum"). 
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Therefore : 

N N-1 

A0 {tN) = ^(fc/3z - k^.).Q iU) - - ^c.^) (23) 

i=l i=l 

which has a very simple form when the origin of coordinates is chosen such that 
qi = 0, and when the Bragg condition pi + = is satisfied. 



8 Conclusion 

In this paper we have used the ABCD^ formulation of atom optics and the 
ttt theorem to establish two theorems valid for a time-dependent Hamiltonian 
at most quadratic in position and momentum operators. The first one gives 
a compact expression of the action difference between two homologous paths. 
The second one gives an analytical expression of the global phase shift for atom 
interferometers in the case of such a Hamiltonian. 

Consequently this analytical expression provides a simple way to calculate 
exactly the phase shift in this case, and then to calculate perturbatively for 
example the effect of a third-order term in the external Hamiltonian (necessary 
for space missions like HYPER 0). For example, one can calculate exactly 
the global phase shift due to gravity plus a gradient of gravity plus a rotation, 
and then calculate perturbatively the effect of a gradient of gradient of gravity. 
These calculations and the application to specific cases (gravimeters, gyrom- 
eters, atomic clocks...) will be detailed in a forthcoming article where we 
recover well-known perturbative results ([^, ||l^) from exact expres- 

sions. 



9 Appendix 1 

In a case of a Hamiltonian at most quadratic in position and momentum op- 
erators, the Hamilton principal functions concerning two pairs of homologous 
points are also at most quadratic in positions (owing to the Hamilton-Jacobi 
equation, see |^): 

Sa {qal,qa2) /ma =a + b.q^l + C.qa2 + qal-d.qal + qal-e.qa2 + qa2-f-qa2 (24) 



Sp (g/31, g/32) /mp = a + h.qpi + c.g/32 + qpi-d.qpi + qi3i.e.qp2 + qp2.f-qp2 (25) 
where a is a scalar, h and c are vectors, and d, e and / are matrices (see 0). 
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We can define Pai, Pa2, P0i, Ppi such that: 
Pal = -Vg„i 

Pa2 = Vq^s 
PHI = -Vg^i 

and obtain the following expression: 

S S \ 1 

— - = 7, {Pa2 +P/32) • {qa2 - 9/32) - (Pal +P0l) ■ (Sal - 9/31) (30) 

The same relation holds for the classical action concerning two actual paths 
with a common drift time (homologous paths). This yields an other demonstra- 
tion of the first theorem expressed in section 4. 
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